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O ' Mean-field theory describes magnetohydrodynamic processes leading to large-scale magnetic fields in various cosmic objects. In this study 
' magnetoconvection and dynamo processes in a rotating spherical shell are considered. Mean fields are defined by azimuthal averaging. 
In the framework of mean-field theory, the coefficients which determine the traditional representation of the mean electromotive force, 
l including derivatives of the mean magnetic field up to the first order, are crucial for analyzing and simulating dynamo action. Two methods 

are developed to extract mean-field coefficients from direct numerical simulations of the mentioned processes. While the first method does 
, not use intrinsic approximations, the second one is based on the second-order correlation approximation. There is satisfying agreement 
T— H ■ of the results of both methods for sufficiently slow fluid motions. Both methods are applied to simulations of rotating magnetoconvection 
^ and a quasi-stationary geodynamo. The mean-fleld induction effects described by these coefficients, e.g. the a-effect, are highly anisotropic 
in both examples. An a-^-mechanism is suggested along with a strong 7-effect operating outside the inner core tangent cylinder. The 
turbulent diffusivity exceeds the molecular one by at least one order of magnitude in the geodynamo example. With the aim to compare 
mean-fleld simulations with corresponding direct numerical simulations, a two-dimensional mean-field model involving all previously 
determined mean-field coefficients was constructed. Various tests with different sets of mean-fleld coefficients reveal their action and 
, significance. In the magnetoconvection and geodynamo examples considered here, the match between direct numerical simulations and 
mean-fleld simulations is only satisfying if a large number of mean-field coefficients are involved. In the magnetoconvection example, the 
azimuthally averaged magnetic field resulting from the numerical simulation is in good agreement with its counterpart in the mean-fleld 
C " ^ model. However, this match is not completely satisfactory in the geodynamo case anymore. Here the traditional representation of the 
mean electromotive force ignoring higher than flrst-order spatial derivatives of the mean magnetic field is no longer a good approximation. 

^ Keywords: Magnetohydrodynamics; Mean-Field Theory; Dynamo Coefficients; Magnetoconvection; Geodynamo 

o ■ 

22 1 Introduction 

i> ' 

.!-H ' Large-scale cosmic magnetic fiel ds as the Ea rth's, the solar and the galactic magnetic fields are maintained 
^ ■ by hydromagnetic dynamos (e.g. WeisdEooii ) . However, global computational dynamo models simulate only 



in 

On 



the geodynamo resonably well, while it remains difficult to tackle the stellar or galactic dynamo problem 
in this way. Despite the increasing computational power, a direct numerical treatment of the governing 
equations is not yet feasible in the latter cases, because of the hu ge range of spatial and temporal scales 
needed to be resolved there ()Tobiasll2n'ni . IWeiss and Tobiaai,2flnft. .Shukurovtl2nr)2i'l . First atte mpts of global 



MHD simulations have been carried out by Oilman dGilman and MilleTJIl 98ll lrTi]ma.nlll98.'j ^ adopting the 
Boussinesq approximation and later by Glatzmaier (|Glatzmaieilll984l . ll98fJ ) using an anelastic model. Both 
aimed at modelling the solar dynamo. Though cyclic dynamo solutions were obtained in some cases, the 
models exh i bited a wrong poleward migration of the magnetic field ( Gl atzmaier 1985 ). In a recent attempt, 
Bnm et al. 1)20041 ) succeeded in simulating a solar-like differential rotation, but could not reproduce basic 



features of the solar cycle yet. 

An a lternative approach is provided by mean-field electrodynamics ()Steenbeck et al ]IT9fi^. lMoffattl lT97^. 



Krause and Radler 1980), which is a theory focussing only on large scale, i.e. averaged fields. Highly 
complex small-scale or residual parts need not to be known in detail, only the averaged cross product of 
the residual velocity and magnetic field, in the following called the mean electromotive force, is relevant 
and accounts for the evolution of the mean field. Advantageously, the difficulties in resolving the small- 
scale structures can be avoided. Usually, the action of the small-scale velocity on the mean magnetic 
field as expressed by the mean electromotive force is parametrised, and the parameters are the mean-field 
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coefficients defined below. Most prominent among them are those which describe the a-effect. They are 
closely related to the fundamental induction effect of cyclonic convection f|Parkeriil955 ). Other mean- field 
coefficients describe the advection and diffusion of the mean magnetic field. In other words, mean-field 
theory supplies theoretical insight as well as formalised physical concepts in order to interpret and, in 
principle, also to quantify dynamo action. 

Despite their rel ative simplicitv. mean-field rnodels have succes sfully reproduced basic features of the 
solar cycle (see e.g. Steenbeck and Krausel 19691 . OssendriiverlbOf),'^ and are nioreover unique in cohe rently 
simulating many features of the magnetic field in spiral galaxies ()Beck et al\ hmd . IShukii7^l2nn2l V But, 
whether or not mean-field models show dynamo action depends strongly on the mean-field coefficients. 
They are not known a priori but instead in general determined in a heuristic way. 

Concerning the geodynamo, the situation is different. Many features of the Earth's magnetic field are 
successfully reproduced by nonlinear three-dimensional simulations of the magnetohydrodynamics in the 
Earth's core. Although some model parameters still do not reach realistic values and in particular viscous 
effects and therefore the size of viscous boundary layers are by orders of magnitude overestimated, the 
simulations exhibit a magnetic field dominat ed by an axial dipole at the Earth's surface that is main- 
tained over several mag netic diffusion times ( Glatzmaier and R,obertslll995^ . iKuang and BloxhamI Il99^ . 
Christensen et al\\l99st \. In addition, the time-dependence of the dipole momen t, including secular vari- 



ations , excursions and reversals, resembles the observed Eart h's magnetic field ( Glatzmaier and Robert^ 
1995bl . iKutzner and ChristensenlEnoi . IWicht and 01sonll2nni ) . The success in simulating the geodynamo 
can be attributed to the rather moderate vigour of turbulence in the Earth's outer core compared to the 
much more turbulent dynamics in the solar convection zone. The difference is formally exp ressed in ter ms 
of very different magnetic Reynolds numbers: While Rm ~ 10^ for the Eart h's outer co re ('Fearn| ll998l ). a 
representative value of Rm at the base of the solar convection zone is 10^° (|OssendriiveL,.2nn3, ). 

Even though present geodynamo models are fully self-consistent, the interpretation of the mechanism 
bv which thev generate magne t ic field relies frequentlv in a heuristic man ner on mean-field concepts 
( Glatzmaier and Roberts! 1995al . Kagevama and Satol 1997 , Olson et al. 19991 ) . Thus mean-field theory is 
very useful and indispensable at the present time. However, the applicability of mean-field concepts as 
tools for analysing dynamo processes in numerical simulations suffers from the poor knowledge of the 
mean-field coefficients and reliable methods to derive them. 

There are two seminal approaches which have been followed in order to determine mean-field coeffi- 
cients. The first one aims at derivin g (quasi-)analvtical expressions of the m ean electromotive force and 
the mean-field coefficients (see e.g. iMoffattI Il97^ . Krause and Radlei 198Ch . This requires the integra- 
tion of the governing equation for the residual magnetic field with the help of closure methods. Most 
commonly used is the second-order correlation approximation (SOCA) in which only statistical moments 
up to the second order are taken into account, while moments of higher order are neglected. Early in- 
vestigations with SOC A and simple assun iptions on the turbulence , starting with the seminal paper by 



Steenbeck et a,l\ \i 



tions on 



, ^ , a re summarize d in ^Ri idleiJ 1 198r , ll99,4Enflflh. More sophisticat e d assu mpti 

the turbulent motion let iKichatinov an d Riid ieeil 1 199/ ) and iRiidiger and Kichatinov (|l993l l provide a- 
coefficie n ts in the high conductivity limit for various rotation rates. In the context of the galactic dynamo, 
Ferrierd (Il99l Il993al lbh considered supernova explosions and derived expressions for the components of 
the a- and /3-tensors in cylindrical geometry. S teps beyond SOCA for s pecific flows were undertaken by 
Nicklaus and Stixl (ll988l 'l. lCarvalhol (Il99lll994) and lSchmittJ (119841 . 1200.11 1. In the context of the Karlsruhe 
dynamo experiment the velocity field possesses a number of symm etries, which allow a rather direct compu- 
tation of mean-field coefficients ( Radler et a/.lll997 . Il998!. '200 2a|l3). A p articular approa ch beyond SOCA 
is the so-called r-approximation, fir s t introduced by IVai nsh tein and Kichatinovl (Il983ll. recently reviv ed 



by Rogachevskii and KleeorinI (l2003h . iRadler et al. J[20o3) and iBrandenburg and SubramanianI |2005 



and critically reviewed bv lRadler and Rheinhardt! (|2006l l 



The second approach to the mean-field coefficients makes use of numerical modelling, in which the mean 
electromotive force, ^, as well as the mean-field, B, are determined numerically as output of a MHD- 
simulation. Further on, a linear relation between S and B is assumed, which has to be inverted in order to 
solve for the unknown mean-field coefficients. A fundamental problem related to this approach is that the 
number of unknown variables is in general much higher than the number of equations resulting from the 
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linear relation between £ and B. Therefore, all work presented so far refers either to specific situations in 
which certain constraints reduce the number of mean- fi eld coefficient s befor ehand, or some of the mean-field 
coefficients are considered as small and are neglected. Zieeler et al. calculated the a-tensor for the 

galactic dynamo on the basis of numerical simulations of supernova explos i ons and con firm ed the results 
bv lFerrier^ ( 1993al ). Following a similar approach, Ossendriiver et all (|2nnil . Enfl3 ) and Kapvla et al. 



()2flnfil ) considered magnetoconvection in the so lar convection zone an d used box simulations to determine 
the a-tensor in Cartesian geometry. Similarly, Giesecke et al. (200^) derived a geodynamo a-effect from 
box simulations of rotating magnetoconvection. 
/3-tensor as well for accretion and gala ctic disks, 



n a further attemp t to determine not only the a.- but the 
Brandenburg and Sokoloff (2002 ) and Kowal et al\ (l200(ih 



applied box simulations of turbulence. ISchrinnen l|2005l ) calculated all mean-field coefficients in a spherical 
shell with a convection pattern relevant for the g eodynamo. The method and results are described in detail 



in this paper. For shorter contributions see also ISchrinner et al. I (|200,1 I200(tI 'I 



The idea of this work is to take advantage of global numerical simulations of rotating magnetoconvection 
and the geodynamo and to compare them with respective mean- field calculations, where mean fields are 
defined by azimuthal averaging. This will lead to an estimation of the reliability of mean-field theory 
and its often used approximations. Furthermore, such a comparison will help to improve the conceptual 
understanding of dynamo mechanisms which are observed in the numerical simulations. 

As already pointed out, both aims are intimately associated with the derivation of the corresponding 
mean-field coefficients. Hence, emphasis is placed on the developement of two methods which contribute 
to each of the principal approches mentioned above. Both methods have been applied to a simulation of 
rotating magnetoconvection and a quasi-stationary geodynamo. They are consistent with each other in a 
parameter regime in which the second-order correlation approximation is justified and serve as powerful 
tools to determine a number of relevant mean-field coefficients. While most of the quoted earlier work 
refers to a Cartesian geometry, global mean-field coefficients for the astrophysically more relevant domain 
of a rotating spherical shell are presented here, and specific problems related to the spherical geometry are 
discussed. 

The plan of the paper is as follows: In section 2 the equations, the boundary conditions and the pa- 
rameters of the considered numerical models are given. Section 3 summarizes the mean-field concept and 
introduces the mean electromotive force. In section 4 two approaches to determine the mean-field coeffi- 
cients are developed, a general numerical one and a semi-analytical approach using SOCA. In section 5 
the mean-field coefficients for two models, rotating magnetoconvection and a quasi-steady geodynamo, are 
derived and their quenching as well as limitations of the validity of SOCA are discussed. A comparison 
of the mean fields for both models, derived from numerical simulations and from mean-field theory, re- 
spectively, is made in section 6, and the range of validity of the representation of the mean electromotive 
force are discussed. Section 7 summarizes our conclusions. In the appendices details of the derivation of 
the mean-field coefficients are given and the mean-field energy balance is considered. 



2 The numerical models considered 

In both models, magnetoconvection and geodynamo, a rotating spherical shell of electrically conducting 
fluid is considered in which the fluid velocity V, the magnetic fleld B and the temperature T are governed 
by the following equations using the Boussinesq approximation: 



+ (V ■ V)V - i/V^V + 2 n X V = --VP + argT + — (V x B) x B 
at Q 

= V X (V X B) + T/V^B 



+ V • VT = kAT 



dT 
'dt 



V-V = V-B = 0. 



(1) 
(2) 

(3) 
(4) 
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Here P means a modified pressure, q is the mass density of tlie fluid and /U its magnetic permeability. 
Further u, rj and k are kinematic viscosity, magnetic diffusivity and thermal conductivity, and ot the 
thermal expansion coefficient, all considered as constants. The motion is measured relative to the uniform 
rotation of the shell with angular velocity d = Qe^ where is the unit vector in the direction of the 
rotation axis. The gravitational acceleration is specified by g = gor/rQ, with being its value at the outer 
boundary r = tq. The ratio of inner to outer radius of the shell is rj/ro = 0.35 and thus the shell width 
D = 0.65 ro for all simulations considered here. The original form of the buoyancy term is axsiT — Tq) 
with To describing the temperature distribution of a reference state. Here Tq is assumed to depend on r 
only. Then a-rgTo can be represented as a gradient, which is absorbed in the pressure term in equation (^. 

For the velocity V no-slip boundary conditions are adopted, V = at r = and r = rp. Moreover, 
all surroundings of the spherical shell are assumed as electrically non-conducting, so that the magnetic 
field B continues as a potential field in both parts exterior to the fluid shell. In the magnetoconvection 
case a toroidal field is imposed via inhomogeneous boundary conditions. The temperature is assumed to 
be constant on the boundaries such that T = 6T at r = rj and T = at r = tq. 

Measuring length, time, magnetic field and temperature in units of D, D'^/u, {gfirii})^/"^ and 6T, the above 
equat ions can be written in a. non-dimensional form which contains only four non-dimensional parameters 
(e.g., IChristensen et a/.ll200lh . These are the Ekman number E, the modified Rayleigh number Ra, the 
Prandtl number Pr, and the magnetic Prandtl number Pm, 



E = u/nD^ , Ra = axgoSTD/vQ. 
Pr = V I K ^ Pm = u/rj . 



(5) 



In order to characterise the results of the simulations, the magnetic Reynolds number Rm and the Elsasser 
number A 



Rm=vD/rj, K = B'^/q^t]VL 



(6) 



with V interpreted as r.m.s. velocity and B as the r.m.s. value of the magnetic field inside the shell are 
used. 

For t he numerical soluti on of the above equations, a code is used which was constructed in its original 
form by Glatzmaiei (|l984h . This version of the code solved the anelas t ic ma gnet qhvdrodvnamic equations 
in a spherical shell to simulate stellar dvnamos. lOlson and Glatzmaiei and lChristensen e,t all 

later applied a modified version of the numerical model to run magnetoconvection and dynamo simulations 
in a rotating spherical shell adopting the Boussinesq approximation. The co de has been validated by 
benchmarking it with other three-dimensional models ( Christensen et al. 200 ll ). 



3 The mean-field concept 

3.1 The mean electromotive force and mean-field coefficients 

Within the scope of this work, the mean-field concept is applied to the induction equation Q only. 
In the following, we refer to a spherical coordinate system (r, i}, ip) whose polar axis coincides with the 
rotation axis of the shell. Mean vector fields are defined by averaging the components of the original fields 
over all values of the azimuthal coordinate (/?, e.g., B = Br{r,'d)er + B^{r,'d)e^ + B^{r,i})e^ in which 
Br{r, B^{r, and B^{r, ??) are the azimuthal averages of Br{r, i}, if), B^{r, i}, (p) and -B,^(r, ??, ip). Note 
that with this definition of mean fields the Reynolds averaging rules apply exactly. Of course, all mean 
fields are axisymmetric about the polar axis. 

Subjecting the induction equation @ to this averaging yields 

a'p 

_ = Vx(VxB) + Vx5 + 7?v2b, (7) 
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with the mean electromotive force 



f: = vxb, (8) 

in which v and b are the residual velocity and magnetic fields, b = B — B and v = V — V. If v and V 
are given, the calculation of S requires the knowledge of b, which is governed by 

f)U 

_ = V X (V X b) + V X (v X B) + V X G + r/V^b, (9) 
at 



with G = vxb — vxb. According to and Q, 5 is a functional of v, V, and B, which is linear in B. 
We adopt the frequently used assumption that b vanishes if B does so. This excludes the possibility of a 
dynamo with B = 0, in other context referred to as "small-scale dynamo". Then £ is not only linear but 
also homogeneous in B and can be expressed in the form 

£,{r,^,t)= [ [ K^x{r,^,t;r','d',t')Bx{r',^',t')dv'dt' (10) 
Jv Jt'<t 

with some kernel Kj^x- Here and in what follows indices like k or A are used for r, -i? or and the summation 
convention is adopted. The integration is over the whole fluid shell, V, and all times t' < t. 

It seems plausible to assume that £ at a given point in space and time depends only on quantities in 
certain surroundings of this point. This implies that K^x differs only for sufficiently small |r — r'\, |i? — 
and t — t' markedly from zero. We adopt here the assumption that B varies only weakly in space and time 
so that its behavior in the relevant surroundings of a given point can be well described by B and its first 
spatial derivatives in this point. This brings us from ()1U|) to 

^■K — CIkX^X + O^Xr^ H OkAi? TT^ ■ i^J-j 

or r on 

Note that all higher than first-order spatial derivatives and all time derivatives of B are ignored. It remains 
to be checked to which extent their neglect is justified for the considered examples. The representation (jllj) 
contains 27 independent coefficients a^A, fe/tAr and 6kAi3) which are determined by v and V and, considered 
as functionals of these quantities, independent of B. We have 

a^x{r,^,t)= [ [ K^x(.r,^,t;r',^',t')dv'dt' . (12) 
Jv Jt'<t 

and analogous relations for b^Ar and fo^Ai? with additional factors (r' — r) or r(i?' — ■(?), respectively, in the 
integrand. 



3.2 A more general representation of the mean electromotive force 

We have introduced the representation (|1H) of the mean electromotive force £ under special conditions. 
In particular we relied on a spherical coordinate system and restricted ourselves to mean magnetic fields 
B which are axisymmetric about the polar axis of this system. 

In general, the mean electromotive force £, if no higher than first-order spatial derivatives and no time 
derivatives of B are taken into account, is written in the form 

5 = aB-hbVB (13) 

with a second-rank tensor a and a third-rank tensor b. This relation is understood as a coordinate- 
independent connection between the vector £, the vector B and its gradient tensor VB. It applies indepen- 
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dently of symmetries of B. In a Cartesian coordinate system ((T^ takes the form Si = SijB j + bij^dB j / dx^- 
When changing to a curvihnear coordinate system dBj/dxk turns into a covariant derivative. 
It is useful to rewrite into the equivalent relation 

5 = -a-B-7xB-/3-(VxB)-5x(VxB)-fv - (VB)^'^^), (14) 



Radlei ( 198Cll . 200(J ) or Radler and Steoanov |2006.) . Here o; and /3 are symmetric second-rank tensors, 



see 

7 and S vectors, k is a third-rank tensor symmetric in the indices connecting it with (VB)^'^^™), the latter 
being the symmetric part of the gradient tensor VB. The relationship between the components of a and 
b and those of a, 7, /3, (5 and k is given below. 

The representation 1)14(1 of £ allows a discussion of the individual induction effects. The a term describes 
the a-effect, which is in general anisotropic. The 7 term corresponds to a transport of mean magnetic flux 
like that by a mean motion of the fluid. The /3 term and also the 6 term can be interpreted by introducing 
an anisotropic mean-field conductivity. The k term covers various other influences on the mean field, which 
are more difficult to interpret. 

In contrast to the 27 coefficients o^A; ^KAr and b^Ai? in ()11|) we have here in general 36 independent 
components of a and b, or of a, (3, 7, d and k. The lower number in the case of (|11() is due to the 
assumed axisymmetry of B. We stress that the d^x, b^xr and ^^Ai? should not be considered as tensor 
components, whereas a and b as well as a, /3, 7, 5 and k, are indeed tensors or vectors with the usual 
defining transformation properties under changes of the coordinate system. 

We may, of course specify the relation H13|) to our spherical coordinate system and to axisymmetric B. 
When doing so in the above sense, that is, taking the covariant forms of the derivatives of B, we arrive at 
relations for of the form ()llj) with 

d^r = Snr + b^m/r + bi^^^/r , d^^ = a^^ - b^ri^/r + cot b^^^/r , 

dKip = - {bnnp + cot 1? bK-ff^) /r , h^Xr = b^Xr , &kAi9 = bi^x-0 ■ (15) 

We may understand (P3|) as a system of 27 equations which determine the a^r + b^p^p/r, a^,? + cot t9 b^^ip/r, 
Snip — {b^rip + cot ^bi^i)^p)/r, b^Xr and b^x^ if the d^x, b^Xr and b^x^ are known. Clearly we have the freedom 
of arbitrarily choosing the b^Xtp- Of course, this choice also influences the a^A- It is, however, without any 
influence on £. For what follows we put simply b^Xtp = 0. Then we have 

a^r = Sfcr bf^^^/r , a^j^ = d/^^ + bf^j.^/r , a^i^ = fl^i^ 

bnXr = bnXr , b^xS = &kAi? , b^X^ = . (16) 

Using this we may also express the components of the a, /3, 7, S and k in the spherical coordinate 
system by the a^^Ai ^^Ar and 6kAi?- The relations of the components of a, 7, /3, 5 and k to that of a and b 
are given by 

"kA = -l/2(aKA + axn) , 7k = l/2eKA/jaA/i 

PkX = '^/^(^KiMubxfiu + eXfiubntiiy) (17) 
= -l/4(faAKA - ^AAk) , l^nXfi = -'^/^bi^Xfi + bi^^x) ■ 

Combining this with (|16() we find 

arr = -{drr - brm/f) , Qr,? = "i^r = "1/2 {dr^ + brr-d/r + a^r " hm/r) 

OLr^p = a^pr = -1/2 (dr^ + d^r " b^i}^/r) , a^^ = -(a^^ + b^r^/r) 



a^ip — ot^p^ — —1/2 {d^^p + d^p^ + 6(^ri?) i « 



ipip (^Lp:p 
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l^rr = —1/2 hrifd , Pr^ = P^r = 1/4 [bnpr — b^ip-&) , Prifi = Ptpr = 1/4 (6rri? — h^p^p{) — hr{)r) 
P-d-d = 1/2 h^^pr , P^^p = P^p^ = 1/4 + b^p^pr — b^^r) i Pipip = 1/2 (bi^ri? — ^</3i?r) 

jr = l/2{di}p - d^i} - b^ri}/r) , 7^ = 1/2 (a,^^ - 6^^^/r - a^^) (18) 
7</3 = 1/2 {ctr^ + brr-d/r - d^r + hm/r) 

6r = 1/4 — bi)r-d + b^^r) , ^i) = 1/4 (6rri9 — br-ffr + ; ^^t/J = — 1/4 (6r(/3r + &i9<^i?) 

^Krr — bj^^^ , I^Kr'd — ^K'&r — 1/2 (^b^ri^ ~1~ bi^{}^^ , l^Krcp — ^Kipr — 1/2 b/^^r 
^K-d'd — ~bK-d-d ) I^k^Ilp — I^KLpi) = —1/2 ^ki^i? ) t^Kipip = . 

4 Determination of mean-field coefficients 

For the determination of the mean- field coefficients o^A; ^^Ar and b^^x^ from the results of the direct 
numerical simulations, two different approaches have been used, which are explained in the following. 

4.1 A numerical approach (approach I) 

We start from equation © for b but interpret B as a steady "test field" Bt- More precisely, we require 
that b satisfies 

— - Vx(Vxb)-VxG - r^V^b = V x (v x Bt) (19) 

inside the conducting shell and continues as a potential field in both parts of its surroundings. The initial 
conditions loose their importance after some transient period. Calculating £ numerically according to (jSJ 
and (|19|) for a given Bt and inserting the result in (jlip provides us with three equations for the wanted 27 
coefficients. We therefore carry out such calculations with the same V and v but nine different test fields 

Bt^ to obtain nine mean electromotive forces i = 1, 2, . . . , 9. As far as the conditions for the validity 
of Hll|) are fulfilled we have then 

^.«=a.A:B« +6«A.^ + 6.A.-^, i = l,2,...,9. (20) 

or r ait 

The (XkA; ^ftAr and ft^Ai? depend on V and v but not on Bt, that is, they do not depend on i. Clearly 
(PU)) represents for each fixed k a set of nine linear equations for the nine coefficients d^x, b^Xr and b^x^. 

The S^'^') and B^x occurring in these equations are then known quantities. Provided the Bt are properly 
chosen, we may solve (1^0]) to find all o^a, ^KAr and 6kAi?- The quantities £ and Bt as well as the d^x, 
b^Xr and b^x^ are in general functions of position and time. The described procedure can be applied to all 
positions and all times. 

In what follows V and v will be extracted from the numerical simulations with the models defined in 
section 121 that is, by the equations Q-Q. Technically, parallel to the numerical solution of a problem as 
defined there, also the nine solutions b of (IT9l) have been calculated. 



There are some constraints on the choice of the test fields Bt . Of course, they have to be axisymmetric. 
They also have to be linearly independent. Otherwise there are no unique solutions of the equations H2U|) . 
We further have to require that all higher than first-order spatial and all time derivatives of the test fields 
are equal to zero, or at least sufficiently close to zero, since otherwise (|2()|) is no longer justified. Fortunately, 

however, unlike B, the test fields B^^ need neither to be solenoidal nor to satisfy any boundary conditions. 
This follows from the fact that relation (|10j) can be derived using nothing else than the definition £ = v x b 
and an equation for b that formally agrees with ©, in which B, however, is understood as any axisymmetric 
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Table 1. A set of test fields B.p used for the determination of d^x a-nd 
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vector field. A set of test fields Brp used in our calculations is given in tabled Note that not all of these 
vector fields are regular at the polar axis. Consequently, V x (v x Bx) could become singular if the axis 
were included in the grid and v were different from zero there. We also experimented with other test fields 
and verified that the mean-field coefficients are independent of their particular choice as long as the above 
constraints are obeyed. 

As explained above, in the determination of the mean electromotive force € often SOCA is used. It is 
defined by cancelling the term with G in ((TO)) . Our procedure for the calculation of the a^A, &ftAr and ft^Ai^ 
also works on this level. 



4.2 A semi- analytical approach using SOCA (approach II) 

We start again from equation @ but introduce some simplifications so that the remaining equation for b 
allows an analytical solution. In that sense we restrict ourselves to the case V = 0. Furthermore we accept 
the second-order correlation approximation and cancel the term with G. Finally we consider only the 
steady case, that is, assume v, b and also B to be independent of time. With these assumptions equation 
© turns into 

ryV^b = -V X (v X B). (21) 
In the solutions of the problems defined in section |2l the velocity v is represented in the form 

V = -V X (r X V0) - r X VV' (22) 

with scalars (j) and given by 

(t> = Y.4>T{r)Yr{^,^), ^l^ = Y.i^T{r)Yr{^.v). (23) 
l,m l,m 

The 0™ and tpj^ are complex functions of r satisfying c/)'^* = and ip"^* = ipf"^, but cf)^ = ipf = 0. 

The Yj^ are spherical harmonics, Yj^{^,(p) = P^''"' (cos ??) exp(im93), with being associated Legendre 
polynomials. In the following the 0™ and ipj^ are considered as given. 

In appendix^the solution of equation l\'21\i for b is derived, that is, b is expressed by the 0™, ipf^, YJ™ 
and the components of B. On this basis £ has been calculated. It occurs at first in a form analogous to 
(fTfl|) . more precisely 

£^{r, [ K,x{r, r' , ^')Bx{r', dv' . (24) 
Jv 

The kernel i^^A is determined by the 0["(r), ^^{r) and the P™(cost?). As an example, K^r is explicitly 
given in appendix 
Knowing the kernel -ftT^A we may calculate the a^A according to 



Jv 



(25) 
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Figure 1. The radial velocity in the magnetoconvection case with 6-fold azimuthal symmetry at r = 0.59 ro, normalized with its 
maximum given by Vr = 16.98 »y/Z). In the grey scale coding, white and black correspond to —1 and +1, that is, downflows and 
upflows, respectively, and the contour lines to ±0.1, ±0.3, ±0.5, ±0.7, ±0.9. 

This relation turns into those for b^Ar or bf^x^ if the factors (r' — r) or r(t?' — t?), respectively, are additionally 
inserted into the integrand. Note that B does not enter the expression for Ki^x and that the mean-field 
coefficients are thus independent of the mean field. 
In fact only the d^x have been calculated so far. It was found, e.g., 



.(r,^)=- y: 



1,1' ; m>0 



/;(r,r')ReUr(r)C(^') 



ri 



/Kr,r')Im(0r(OC(^') - ~gi{r,r')Mi>r {rW{r 



r dr 



xm QITJ(^)+QIIi(^) /sin^ , 



(26) 



where /;(r, r') and gi{r, r') are Green's functions, and Rf^i{'d) and Q",]{^) specific combinations of associated 
Legendre polynomials P[^{cos'd), all explained in appendix 1X1 The other a^A ^.re also given there. 



5 Mean-field coefficients for rotating magnetoconvection and a geodynamo model 
5.1 Rotating magnetoconvection 

The first example considered is adopted from lOlson et al. with the governing parameters E 



3 X 10 ^, Ra = 94(= l.hRac, where Rac means the critical value of Ra) and Pr = Pm = 1. Moreover, an 
axisymmetric toroidal magnetic field is imposed via an inhomogeneous boundary condition of the form 

= Bo sin(2t?) at r = ri, r = tq . (27) 

Figure n shows the radial velocity field at r = 0.59 tq for an Elsasser number of the imposed field Aq = 1, 
where Aq is defined according to (6) but with B replaced by i^o- A typical columnar convection pattern is 
revealed. Apart from a steady azimuthal drift of the convection columns the fiow and also the magnetic 
field are steady. In addition, the azimuthal drift and so the remaining time dependence can be removed 
by a transformation to a corotating fra me of reference. The electro motive force £ is identical in both the 
original and the corotating frame (e.g., Riidler and Stepanov 200()h . 



The magnetic Reynolds number Rm is about 12, that is too low for the onset of self-sustaining dynamo 
action. Nevertheless, fundamental effects of the convection, which are of high interest for a dynamo, can 
be analysed in terms of mean fields, for example the generation of a poloidal from a toroidal field. The 
ratio of magnetic to kinetic energy density is about 10, and the resulting field strength is described by 
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Figure 2. Components of the symmetric cc-tensor and the -('-vector in a meridional plane in the magnetoconvection example, 
determined by method I, in units of r]/D. For each component the grey scale (white - negative, black - positive values) is separately 
adjusted to its maximum modulus. Note the negative signs in the definitions of ol and -y in equation 1141 . 

A ^ 0.6. 

Figure 121 shows the six independent components of the tensor a and the three components of the vector 
7 in a meridional plane derived by approach I. All mean- field coefficients are essentially determined by 
the columnar convection outside the inner core tangent cylinder. As a consequence of the symmetry 
properties of the velocity field and the induction equation, all mean-field coefficients are either symmetric 
or antisymmetric with respect to the equatorial plane. The diagonal components of a, for instance, are 
antisymmetric, in its major contributions negative in the northern and positive in the southern hemisphere. 
Among the Q-components, a^pip dominates, indicating that the generation of a poloidal field from a toroidal 
one is more effective than the reversed process. However, due to the other non-vanishing components, 
especially Urr, Ori? and a^^, generation of toroidal field by an a-effect also takes place. 

The mean-field diffusivity tensor D is given by 



(28) 



Its components are shown in figurelHl Although the molecular magnetic diffusivity is rather high (Pm = 1) 
and the vigour of the convection is rather low {Ra = 1.5Rac), the turbulent diffusion is of the same order 
as the molecular one in the convection region. 

The diffusivity tensor D has the interesting property of being positive definite everywhere. This can 
be concluded from the fact that all its diagonal elements as well as all sub-determinants are positive. As 
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Figure 3. The symmetric diffusivity tensor D in the meridional plane in the magnetoconvection example, determined by method I, in 

units of T). Grey scales as in figure 2. 
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Figure 4. (Q!ip(^)rms in units of rj/D (left), the ratio {«r45)rms/{oip¥))rms (middle), and the diagonal components of (/3KA)rms in units of 
r] (right) in the magnetoconvection case, in dependence on the Elsasser number A. 

explained in appendix^ this property implies that the induction effects expressed by D do not contribute 
to a growth of the total magnetic energy stored in the mean magnetic field but favour its dissipation. 
Because the fact that P^x has been defined with some arbitrariness, this statement has, however, to be 
considered with some caution (see also appendix IB)) . 

The d-vector and the ft-tensor (which are not displayed here) have been derived as well, and they are 
used in the magnetoconvection and dynamo calculations of section |H1 



5.2 Quenching of mean-field coefficients 

We point out that the velocity v needed for the determination of the mean- field coeffi-cicnts fl^Ai ^kXt 

and 

6kAi? were taken from simulations with non-zero mean magnetic field B. Therefore, the resulting coefficients 
are already subject to a magnetic quenching corresponding to this magnetic field. In this respect, results 
obtained with various B, measured by the Elsasser number A, are of interest. FigureE](left) shows {a^^)rms, 
the r.m.s. value of a^p^, as a function of A. The increase of this quantity in the presence of a weak magnetic 
field, that i s for small A . is due to an increasing vigor of convection by the relaxation of the geostrophic 
constraint (|Fearnlll998h . A strong magnetic field however inhibits convection and reduces (tti^y::))rms* 

The 

kinetic energy of the convection varies with A similar to {a^p)rms- The other a components are also 
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Figure 5. Components of the symmetric ct-tensor and the •7-vector in a meridional plane in the geodynamo example, determined by 

method I, in units of ri/D. Grey scales as in figure 2. 

quenched. The quenching is however not the same for different components, leading to varying amphtude 
relations among these components for varying strength of the mean magnetic field (figure |1] middle) . In 
addition to the a-quenching, e.g., also a /3-quenching takes place (figure fright). 



5.3 A quasi-steady geodynamo 

As a further example a quasi- steadv geodyna mo model is examined, which has been used before as a 
numerical dynamo benchmark ([Christensen et al. 200i|). The governing parameters have been chosen to 



be £" = 10~^, Ra = 100 (= 1.79i?ac), Pr = 1 and Pm = 5. The columnar convection pattern is similar 
to that in the magnetoconvection example (figure but with a natural 4- fold azimuthal symmetry. The 
intensity of the fluid motion is characterised by Rm ~ 40, and the magnetic energy density exceeds the 
kinetic one by a factor of 20. Again, except for an azimuthal drift of the convection columns, the velocity 
field is stationary. 

The components of the ct-tensor and the '/-vector are shown in figure El Among the ct-components, 
a^pip again dominates, indicating a very efficient generation of poloidal from toroidal magnetic field. The 
components Q^r, c^ri? and are somewhat lower in amplitude. Since the infiuence of the differential 
rotation on the generation of toroidal field is negligible, this example can be classified as an a^-dynamo. 
The imbalance in the amplitudes of the a-components is refiected in the larger strength of the mean 
poloidal field compared to the mean toroidal field. As before in the example of magnetoconvection, the 
7-effect acts to expel flux from the central dynamo region where the convection takes place. 
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Figure 6. Components of the symmetric difTusivity tensor D in a meridional plane in the geodynamo example, determined by method 

I, in units of rj. Grey scales as in figure 2. 
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Figure 7. Components of the 5-vector in a meridional plane in the geodynamo example, determined by method I, in units of rj. Grey 

scales as in figure 2. 



The corresponding diffusivity tensor D is displayed in figure IHl The turbulent diffusivity /3 exceeds the 
molecular one, rj, by more than a factor of 10 in the convection region outside the inner core tangent 
cylinder, leading to a very efficient diffusion of the mean magnetic field. 

There is a weak negative contribution to D^p^ at the inner boundary close to the equator, which is 
negligible in the mean- field model of section HOI The diffusivity tensor D thus slightly deviates from being 
positive definite. We argue in appendix El that another than our particular choice 6kA</3 = will lead 
to another D without changing S and thus the physical situation. Furthermore, a parametrisation of € 
considering higher than first-order derivatives of B will also lead to changes of the low-order mean-field 
coefficients. In section we will see the need of a better parametrisation of £ in the geodynamo case. 

As for the (5-effect we recall th a ,t the combination of this effect with a mea n rotational shear may 
constitute a dynamo dRadlerlll Qfij ll 97nl . FmS^ . lEobert jll 971 Ifiadler f^rfi.l l2nO.'^). A dynamo mechanism 
of that kind, however, can not play a dominant role here, since neither V nor 7 imply a sufficiently strong 
shear. As we will see in section 16.51 below, however, the S terms, which are displayed in figure may 
diminish the decay of a mean magnetic field. 
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Figure 8. The quantity {a^ip)ruis in the geodynamo case, in units of rj/ D, with a^p^ determined by methods I and II (solid and dashed 

line, respectively) in dependence on Rm. 

5.4 Limitation of SOCA 

The mean-field coefficients determined by approach I and approach II (SOCA) show for all Rm an almost 
perfect congruence of their profiles. However, mean-field coefficients determined by means of SOCA exhibit 
typically overestimated amplitudes for Rm > 10. 

In figure IHl (ai^i/3)rms is plotted versus Rm. In all calculations from which the mean- field coefficients 
were derived, v was the same, and the variation of Rm is only due to a variation of r] or Pm which, for 
the calculation of b in equations (|19|1 or (|2H) . can be chosen differently from their values in the original 
model simulations. For small Rm the results of approaches I and II coincide and vary linearly with Rm. 
For Rm > 10, the slope of {a^^)rms determined by approach I flattens. In particular, derived by 
approach I, that is, without restriction to SOCA, leads to amplitudes which are, e.g., for Rm = 40 about 
30% smaller than those gained by approach II, that is, with SOCA calculations. The consequences for the 
dynamo action in a mean-field model are studied in the following section. 

Let us add some explanation for the linear dependence of (a(^(^)rms on Rm observed in both approaches 
for small Rm, and in approach II for all Rm. In SOCA, when assuming a steady velocity v, an ct-component 
like Oipip, say simply a, is given by a = fv'^D/rj, where / is a purely numerical factor. We may write this 
also in the forms a = {fr]/D)Rm'^ or a = fvRm. If rj and D are fixed, a appears to be proportional 
to i?m^, which may then vary with v. If, however, v is fixed, a proves to be proportional to Rm, which 
may vary with rj. This corresponds to the results presented in figure |S1 The deviation of the results for 
(a,^(^)rms obtained in approach I from the linearity in Rm indicates that we are no longer in the range of 
applicability of SOCA or that the time variation of v is no longer sufficiently weak. 

The usually given sufficient condition for the applicability of SOCA in the limit of steady motions reads 
Rm <C 1, where Rm is defined with a typical length of the fluid flow. Even if this length is slightly 
overestimated by D used in our deflnition of Rm, our flnding of the applicability of SOCA for Rm < 10 is 
very remarkable. 

6 Comparison between numerical simulations and mean-field models 
6.1 Mean-field model 

In order to compare direct numerical simulations and mean-fleld theory, an axisymmetric mean-fleld dy- 
namo model involving all 27 iiiGa,n- field, coefficients g-kA? ^kXt: 

and 6kAi? has been constructed. The model 

also enables isolating certain dynamo processes. 
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Decomposing the axisymmetric mean magnetic field B in its poloidal and toroidal parts, 

B = Bpoi + Btor (29) 

with 

Bpoi = V X ^e^ , Btor = Be^, (30) 

where e<^ is the unit vector in azimuthal direction, we may write the induction equation ((T)) as 
gj[ I 

Vpoi • V(r sin ^A) + £^ - 7]MA (31) 



dt r sin -d 

dB . „— „f B \ 1 9(F<n/rsini?, rsini?^) ,„ ^ , 



where A' = A — l/(r sint?)^. Here, the notations £poi = {Eri^'d-,^) and Vpoi = {Vr, V^,Q) have been used. 
£ in its dependence on B has to be taken from The above equations are then solved in a spherical 
shell with electrically insulating inner and outer surroundings. Thus, the mean magnetic field is assumed 
to continue as a potential field in both parts exterior to the fluid shell. 

The two coupled equations (|31 |) -(|32 |) are solved by a finite difference method on an equidistant grid 
in radial and latitudinal direc tion. An alternating direction implicit scheme for parabolic equations with 
mixed derivatives according to iMcKee et al. has been used to discretise the equations. This enables 

an efficient implicit treatment of advection and diffusion terms, whereas mixed and higher order derivatives 
are treated explicitly. 



6.2 Magnetoconvection 

How well do the results given by mean-field models match with the corresponding azimuthally averaged 
fields determined by direct numerical simulations? Let us first consider the rotating magnetoconvection 
model discussed in section I^TTl Figure 1^1 presents a comparison between direct numerical simulations and 
mean-field calculations. In the first row, the azimuthally averaged magnetic field components resulting 
from the direct numerical simulation are shown. They correspond in great detail to the results of our 
mean-field model (second row), in which all 27 mean-field coefficients have been used. The poloidal field is 
dipolar with inverse flux spots near the equatorial plane, and the applied azimuthal field is expelled from 
the region occupied by the convection columns. 

A mean-field simulation relying on mean-field coefficients derived in SOCA (third row in figure IHI) fits 
equally well. This reflects that mean-field coefficients given by SOCA, even overestimated by a few per 
cent in their amplitudes though, still lead to a reliable parametrisation of the mean electromotive force 
in this parameter regime. Moreover, amplitude deviations simultaneous in a and /3 might not strongly 
influence the efficiency of the generation of poloidal from toroidal magnetic field and vice versa, as suggested 
by a simple scaling analysis: The efficiency of these processes can be expressed by the dimensionless 
number P = D"^ / I3q. Here, oq and /3o mean typical values for the a-effect and the turbulent diffusivity, 
respectively, and D stands for a typical length scale. Since cx and (3 are likewise overestimated in their 
amplitudes, this factor cancels out and has no influence on P. As a consequence, the resulting field resembles 
the mean field displayed in the second row, even though the applied mean-field coefficients have larger 
amplitudes. 

The mean field components shown in the last row of figurelUhave been determined applying the isotropic 
approximation = o/^ak, /3ak = Pi^Xk- Concerning the co efficients aj and 0t we rely on SOCA results 
for isotropic turbulence in the limit of steady motion (see, e.g.. lKrause and Radlei ( 19801 ) or Radler (|2000t l') 



February 5, 2008 13:24 Geophysical and Astrophysical Fluid Dynamics schrinner 



16 



M. Schrinner, K.-H. Rddler, D. Schmitt, M. Rheinhardt and U. R. Christensen 








itaifjo 








Ufc DM 





Figure 9. Comparison between imiiierical simulations and mean-field calculations in the example of magnetoconvection: azimuthally 
averaged magnetic field components resulting from a direct numerical sinmlation (first row), results given by the mean- field model 
based on all 27 coefficients derived by approach I (second row), mean-field calculation with coefficients derived applying SOCA 
(approach I with G = 0, third row), and mean-field calculation with coefficients for isotropic turbulence (last row). Maxima and 
minima of the field components are given in units of (gurifl)^^^ . Grey scales as in figure 2. 

and put 



ai = -—a - (V X a) 

Sri 



and Pi = 



3rj 



(33) 



where a is the vector potential of v specified by V • a = 0. With this choice of the mean-field coefficients 
the profile of the toroidal field clearly deviates from that in the cases considered before. It is less diffused 
at midlatitudes and mid radii where convection takes place. This difference can be attributed to the 
absence of the 7-efFect. Already in this simple example the isotropic approximation fails to reproduce the 
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Figure 10. Comparison between numerical simulations and mean-field calculations in the example of the geodynamo: azimuthally 
averaged magnetic field components resulting from a direct numerical simulation (first row), results as given by mean- field modelling 
with coefficients derived by approach I (second row), mean-field calculation with coefficients in SOCA (approach I with G = 0, third 
row), and with coefficients for isotropic turbulence (last row). Different from the solution of the direct numerical calculation which is 

stationary, all moan-fiold solutions decay exponentially with decay rates A = 3.5r]/D^ (second row), A = 6.5r]/D^ (third row), and 

A ^ 75r]/D^ (last row). 

axisymmetric field in satisfactory agreement with corresponding numerical simulations. This indicates that 
in general more mean-field coefficients must be taken into account in order to grasp all relevant dynamo 
effects. In addition there are deviations of about 50% in the amplitudes of the poloidal field. 

6.3 Geodynamo 

Consider now again the geodynamo model of section 15.31 In figure ^1 the azimuthally averaged field 
components resulting from the numerical simulation are shown in comparison with results given by mean- 
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, second row: 



£}'^^^ , £^^^ , £^^^ . Maxima and minima are given in units of {r}/D) (g^rjQ)^^^ . 



field modelling. Figure lTUl is organised in the same way as figurelHlbefore. That is, azimuthally averaged field 
components resulting from a direct numerical simulation have been plotted in the first row, the second row 
shows results obtained by corresponding mean-field calculations, the third row contributes results obtained 
by mean-field modelling with the coefficients determined in SOCA, while for the results presented in the 
last row, the isotropic approximation H33|) has been applied. Note that only the direct numerical simulation 
results in a steady dynamo. All mean-field models shown in comparison are subcritical and the magnetic 
fields decay according to B = Boexp(— At), where Bq denotes the field configuration reached after an 
initial transition phase, and the decay rate A is positive in these examples. Therefore, decay rates rather 
than amplitudes are compared. 

As in the previous example, both mean-field models relying on all 27 mean-field coefficients (second and 
third row in figure UJ^ correspond best to the direct numerical simulation and succeed in reproducing all 
essential features of the field given in the first row. However, both mean-field models are slightly subcritical 
with A = 3.5r]/D^ and A = 6.57]/ D^, respectively, which comes along with topological differences in Bip. 
The flux bundles at low latitudes near the outer boundary are more strongly diffused in the mean-field 
models. The high diffusion in this region is due to the strong 7-effect, which leads to an advection of 
oppositely oriented mean toroidal fields towards the equator, resulting in large gradients. 

Although SOCA is, strictly speaking, not justified anymore, the resulting mean field components in 
the third row are remarkably similar to those obtained by mean-field modelling without applying SOCA 
(second row). For an explanation we refer again to the scaling argument given above in the context of the 
magnetoconvection example. 

Again, mean-field coefficients in the isotropic approximation do not lead to reliable results anymore, as 
can be seen from the last row in figure ITUl There are not only differences in the field distribution, but also 
the decay rate, A ^ 75r]/D'^, is drastically high. 



6.4 Limits of the representation of the mean electromotive force 

The difficulties of mean-field models in accurately reproducing mean magnetic fields compared to direct 
numerical simulations, which arise in the example of the geodynamo model, are due to an inadequate 
parametrisation of E. Let us compare f^^'^^, that is £ as immediately extracted from the direct numerical 
simulation, with 5^^^ defined by £"^^1 = d^xBx + b^xr dBx/dr + l/rbf^x-d dBx/d'd. For the computation 
of £^^^, both B and its gradient have also been taken from the direct numerical simulation. We further 
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define the quantity 



< |£DNS_^MF1| > 



^£ = ' -,.DNSi- ' (34) 



where < • • • > means spatial averaging. 

Consider first the magnetoconvection example. Figure HU shows that 5°^^ and S^^^ are in reasonable 
agreement. We find 5£ ~ 0.28. Therefore, a parametrisation of £ considering no higher than first-order 
derivatives of B is adequate in this example. 

In contrast, this is no longer true for the geodynamo model. In this case 5£ > 4 has been found, indicating 
that a parametrisation according to no longer describes the actual € reasonably well. The assumption 
of a sufficiently weak spatial variation of B, which is needed to truncate the series expansion of £ in (jllj) . 
breaks down. In fact, higher order derivatives of B become large and spoil the representation (|1H) of £ in 
a rather uncontrolled manner. 

This finding is consistent with the following observations. The power spectrum of the radially averaged 
mean magnetic field in the example of the geodynamo model exhibits a peak at Z = 4 containing 15% of the 
total power. In contrast, the corresponding spectrum in the magnetoconvection example does not possess 
noticeable contributions for I > 2, suggesting that the mean-field is indeed smoother in this example. 
Furthermore, the magnetoconvection simulation has been repeated with more complicated imposed toroidal 
fields in order to trigger steeper gradients in the mean field. In this way, it is indeed possible to destroy 
the close match between g™^ and £^^^ as seen in figure ITTl 



6.5 Significance of mean-field coefficients 

In order to investigate the significance of the various mean-field coefficients we carried out a number of 
test calculations with different sets of mean-field coefficients for the geodynamo case. As already presented 
in section 6.3, the inclusion of all a^A and bf^xu leads to the decay rate of A = S-brj/D"^. 

In one series of calculations we only used the a^A and disregarded the b^xu, increased however the 
molecular diffusivity to Pm = 1. With all o^a the decay rate of the dominant dipolar mode is A = 2.1 rj/D"^ . 
Using only the diagonal a^K results in A = —A.Srj/D'^. Besides the diagonal terms, drip and a^^ are most 
important. They provide a strong 7-effect, which as already stated above, acts to expel flux from the 
central dynamo region and results in a decay with A = 4.9?7/D^. As seen in section IH^ bv comparing with 
the numerical simulations, the 7-effect is crucial for the geodynamo. Furthermore it leads to a preference 
of dipolar modes compared to quadrupolar modes, that is, modes being antisymmetric or symmetric, 
respectively, about the equatorial plane. 

With all d^x coefficients included, we next studied the influence of the b^xu coefficients, now again 
for a molecular diffusivity described by Pm = 5. With the diagonal terms only, diffusion is significantly 
enhanced and the resulting dynamo solution is markedly subcritical with A = Ibrj/D^. However, not all 
components of b^xu Sire conducive to the turbulent diffusion. If in addition all coefficients which contribute 
to the (5-effect are considered the decay rate decreases to A = 3.0 r]/D^ and the solution gets close to the 
case where all coefficients are included. This might be due to a constructive action of the combination of 
J-effect and some kind of mean shear occurring with V or 7, which, if stronger, could lead to a dynamo 
even in the absence of the a-effect. As explained in appendixEl the <5-effect has no direct influence on the 
energy stored in the mean magnetic fleld. 

Altogether, the test calculations suggest a large set of mean-field coefficients to be considered in order 
to achieve reasonable agreement between mean-field models and numerical simulations. These are the o: 
and 7 terms as well as the (3 and S terms. The k terms, on the other hand, seem not to be of importance. 
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7 Conclusions 

The knowledge of the mean-field coefficients is decisive in order to analyse and to model dynamo action in 
many astrophysical bodies. In this paper, two approaches to determine mean-field coefficients have been 
developed. While the numerical approach does not use intrinsic approximations, the analytical approach 
is based on the second-order correlation approximation. 

The mean-field view is applied to two examples: a simulation of rotating magnetoconvection and of 
a quasi-stationary geodynamo. In both examples similar processes take place: the mutual generation of 
poloidal and toroidal magnetic fields by an a-effect, fiux expulsion from the dynamo region due to a 
7-effect, and a strong turbulent diffusion, which might be moderated by a (^-effect. 

The calculation of mean-field coefficients provides insight into the reliability of frequently applied ap- 
proximations in the framework of mean-field theory. Most dubious among them is the reduction of the 
a-tensor to an isotropic tensor, which leaves dominating non-diagonal components unconsidered. A further, 
important simplification is the second-order correlation approximation. It typically leads to overestimated 
amplitudes of mean-field coefficients, whereas their profiles are rather unaffected. 

The mean-field picture of geodynamo models is completed by the simulation of axisymmetric fields 
by means of a mean-field model, involving all mean-field coefficients determined. Test calculations with 
different sets of mean-field coefficients confirm their relation to the above mentioned dynamo processes. 
In addition, a comparison with azimuthally averaged fields resulting from direct numerical simulations 
reveals their significance. In the magnetoconvection and the geodynamo example considered here, the 
match between direct numerical simulations and mean-field simulations is good only if a large number of 
mean-field coefficients are involved which contribute to a, 7, /3 and S. The application of corresponding 
mean-field coefficients derived in the second-order correlation approximation leads to similar results. 

The reliability of mean-field models relies on a proper parametrisation of £ in terms of the mean magnetic 
field. In the magnetoconvection example the traditional representation of the mean electromotive force 
considering no higher than first order derivatives of the mean magnetic field is valid. However, already 
in the geodynamo example it seems no longer to be justified and the assumption of scale separation is 
not fulfilled. This limits the applicability of some approximations commonly used within the mean-field 
theory. Nonetheless, even in the geodynamo example the spatial structure of the axisymmetric fields 
obtained by mean-field modelling corresponds roughly with that of the azimuthal averages extracted from 
the corresponding direct numerical simulation. 



Appendix A: Approach II derivations and further results 

We first determine b so that it satisfies equation (|21|) inside the fiuid shell and continues as a potential 
field in its inner and outer surroundings. We represent b in the same form as v in (|22|) by writing 



where and Tp{r) are complex functions satisfying a^* = crf"^ and r™* = Tj~^, but cj' = rj^ = 0, 

and Yi^('d,(p) stands for spherical harmonics as explained above. 
We note that ()A1|) implies 




(Al) 



and expanding the scalars a and r according to 




(A2) 



LV = -r-b, LV = -r-(Vxb), 



(A3) 
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-2. 15 f^:.^^adf\ , 1 d^f 



where 

sm if ov \ ov J sva^ -d Oip'^ 

Of course, if b in (|A1|) is replaced, e.g., by V^b, ()A3|) applies with the same replacement. We further note 
that the Y^^ satisfy the eigenvalue equation 

L^Vr + l{l + l)Yr = (A5) 

and the orthogonality relation 

j Yn^, V.) V^) dn = ^0L±}^^5a'6rnn^' , (A6) 

where dfi means sin.'d di!) d(p and the integration is over all "i? and (/? of the full solid angle. Here Ferrer's 
normalization of the associated Legendre polynomials has been adopted. 
Using these relations equation (PT|) can be reduced to 

l^,{rar)-^-^ar = Fr, l^,{rrn-^-^rr = Gr, (A7) 
applying in the shell ri < r < tq, and 

= /^!,t w7 '"i^'^'m /(v X B) • (r X VYDdn 
' Anrjlil + + \m\)\ J ^ ^ ^ ' 

{2l + l){l-\m\)\ 
ATTr]l{l + !)(/ + |m|)! 

again with integrations over the full solid angle. The continuation of b as a potential field in the regions 
inside and outside the conducting shell requires 

— ^ ar = Tr = at r = ri, + = ^'^ = at r = tq . (A9) 

dr r dr r 

The solutions of HA7|) satisfying ()A9|1 can be written in the form 

f-ro pro 

<yT = - fi{r,r')Fr{ry'dr' , = - / r')Gr(r')r' V (AlO) 

Jvi Jri 

with Green's functions // and gi defined by 

— I ^ r' < r 

flir,r')= <; '''^ (All) 



= , '"T^^'u, / (V X (v X B)) • (r X VYD dQ , (A8) 



{21 + l)r 



l+i 

, r' > r 
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and 



gi{r,r') = < 



21+V 



7) h{ny) + -i^h{ny) 



r' r 



-) 

roj 



21+V 



-1 



(A12) 



Mr,r')-(-] Mro,r') 



^0 



r' > r . 



Indeed it becomes clear by inserting of HA10|) that ()A7p is satisfied. Since ()A11|1 and HA12|) imply 

^^^^^^^ - -fl{r,r') = gi{r,r') = at r = n 
or r 

Z \ Jl{r,r ) = gi{r,r ) = at r = tq , 

or r 



(A13) 



()Ain|) also satisfies 

Inserting now the representation of v as given by (|22jl and (|23|) into the relation for given by 
and then the result for into the expression for o"™ in (|A10|) . we arrive at 



<(r) 



(2/ + l)(/- |m|)! 
47r7?/(/ + 1)(Z + \m\y. 



Y^J Mr,r')^[4>P{r')ir^{^') 
-im^lf(r')(gF}(^?') +QI?W)/sin^']s,(r',^') 



r 



+im^^^^-tll</.[r(r') (Ppii^')/ sin B^[r' , \dv' . 



(A14) 



The integration is over the whole fluid shell. When proceeding analogously with the relation for G™ in 
()A8|) . the result for contains derivatives of Br, B^, and B^ with respect to r and -i?. We may remove 
them by means of integration by parts. In this way we find 



rrir) 



(2/ + l)(/- |m|)! 
ATrr]l{l + + |m|) 



1 



+iml'{l' + l)gi{r,r') (l)f{r')PlJi{^')/ sin^' B^{r',^ 
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(A15) 



with gi{r,r') = (l/r') d{r'gi(r,r'))/dr' . 

The relations (|A1|) and HA2|) together with HA14|) and (|A15|) represent the wanted solution of the equation 
(EU) for b. 

Let us now proceed to 5 = v x b. Expressing v according to and (|23() and b according to (|A1|) and 
M we find 



■2 E 



; m>0 



/'(/' + l)Re( ) Qn - Kl + l)Re( V^r^r Qw 



1,1' ; m>0 

+m (lii + i)imUr<^r ) + ni' + i)im( ^r-^r ) i ppi/sm^ 



(A16) 



9^ 



1,1' ; m>0 

The CT/" and tJ^ depend of course on r. The Ppi, and RP'i are functions of i? defined by 



PlTi = pl,'"l(cosi?)pl"l(cosi?) 



dpl"'l(cos^) 
d^9 



(A17) 



^I'l 



dP, (cOS^)dP, (COS^) ^ . 

sin'^T? ' ' i V 



dt? 



We may now insert the results (|A14|) and (|A15|) for and into (|A16|) . Then the take indeed the 
form (|24j) . Unfortunately the iT^A ai'e rather complex expressions. As an example we mention 



K„(r,i9;r',i?') 



(2/ + l)(/- |m|)! 



— y 



/,(r,r')Im(^r(0^J^(O)+5/(r,r')Im(C*(r)C(r') 



+ 



/Kr,r')Re 0rWV'I^(O -5Kr-,r')Re V'r(OC(r-' 
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(A18) 



where fi{r,r') = {1/r) d{rfi{r,r'))/dr. 
Using now H25|) and the orthogonahty relations 



(2/ + - |m|)! 

Dm/ ON • a. a 2l{l + + \m\)l 

(2/ + - |m|)! 



(A19) 



we find drr as given by ((^ . 

In the same way all other dj^x can be calculated. Unfortunately in the cases of a«;^ and a^^ the integration 
over 1? can not be carried out with taking benefit of orthogonality relations like (|A19|) . and the need of 
numerical integrations remains. The results read 



d^rir, ■!?) 



2 



.,.';m>0 I 

- l{l + 1) j^J fi{r, rORe(^r(OC(^')) r'^dr' QJ?(.?) 
+ m^'" /'(/' + l)/Kr,r')Im(</.r(r)C(0' 
+/(/ + l)fi{r,r')lm(^'P*{r)4>r{r')) r'^dr' P;7?(t9)/ sin 



(A20) 



1,1' ]m>0 I 



ri 



+ m 



/'(/' + l)ffz(r,rOlm(c*(OC(r') 
/(/ + l)/i(r,r')Imte*(?')C(^'))l ^'^dr' / sin 



(A21) 
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n (rS)-^ V- {2l + m-\m\)\ 



m 2 \^l{l + l)9i{ry)lm(^4>r{rW'{r')) 
+ l)gi{r, r')lm(0r(r)0^(r')) | r'dr' £ Pll^ii'&'W 
-l{l + 1) J\i{r,r')Re($];*{r)^]li{r')') r'dr' £ Q^„{^')sm^'d^' 

+ 1"" fi{r,r')Re(^^rirWir')) r'dr' f Q^,i{^')sin^'d^' 



+m(Qf0) + Q^,{^))/sm^ m£ |^(/ + l)9i{r,r')Re{4^r{rW'ir')) 

+ iygi{r. r')Re(^r(OC(^')) } r'dr' £ Pl^0)d^' 
+ 1) £ gl{r,r')lm(^l;r{r)^'P{r')) r'dr' £ Q^„{'d') sin^'d^' 
+ 1) £ fi{r, r')Im(c*(0</'P(^')) r'dr' 
X Q'^,i{^')sm^'d'd' \ 



(A22) 



~ , 1 {2l + m-\m\)\ 



x<^ i'{i' + i)Q];i{^) 



KZ + l)5/(r,r')Im(0rWC(r')) 
+ + l)5Kr,rOlm(^r(0</'r(r'')) } r'dr' £ P]J}i{'d')d^' 
+ l{l + l)£ gi (r, rORe (c* (r ) VP (r')) r'dr' £ QJ?, (^') sin 7?'d^' 



^smd'd-d' 



£ 1^1(1 + l)Mr,r')lm(^^rir)4>T'ir') 



Jo 

+ l'{l' + l)/;(r,r')Im((/.p*(r)(/.p(r')) | r'dr' Ppj {'d) / sin 
- l{l + 1) r /Kr,r')Re(c*(r)C(r')) r'dr'QI?(^?) 



(A23) 
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(2/ + l)(/-|m|)! 



x|mZ'(/' + l)P,7?(??)/sint? 



l{l + l)9l{ry)Re(^ct>r{rW'{r')) 

"{I" + l)gi{ry)Re[cl>rirW{r')) } r'dr' £ PjPii'&'W 

- l{l + 1) £ gi{r, r')Im(^r(OV'P(^')) r'dr' £ QI?,(,?') sin^'d^' 

+ /"(/" + 1) r QPK'^') sini?'di?' 
Jo 

l{l + l)Ppi{'&)/s\n^ £ /Kr,r')lin(vr(0'^p(^')) r'dr' 
l'{l' + l)Q^i{^) £ Mr, r')Re(4>r{r)4>7^{r')) r'dr' 
+ l{l + 1)QIJ^ (^) Mr, r')ReUr{rWir')) r'dr' 

Jri 



(A24) 



ar^{r,"9) = -- ^ 



(2^ + l)(;- |m|)! 



isint?'dt?' 



m£ |/(/ + l)gi{ry)lu,[^r{r)^l^{r')) 

- I" {I" + l)/,(r, r')Im(vr(^)'^r(^')) } r'dr' £ PIJ^i{'&')d-d' 
+ l{l + 1) £ gi{r, r')Re{4>r{rW'ir')) r'dr' £ QJ?,(t?') : 

- + 1) £ ~gi{r,r')Re[4>rir)^l^{r')) r'dr' £ QPi{^')sm^'d^' 
m(Q];i{^) + Q^,{^))/sm^ m£ + l)giir,r')Re(^rir)^l^ir')) 
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^"^^,v^.>o^(^ + i)(^ + H)' 
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Appendix B: Mean- field energy balance 

Above we have applied the mean- field concept to the induction equation ((2)). For the following consider- 
ations it is more convenient to start from Maxwell's equations in the quasi-steady approximation. Their 
mean-field version reads 

VxE = -— , VxB = /ij, V-B = 0. (Bl) 

Consider a finite fluid body embedded in free space and assume that there are no causes of B at infinity. 
Then, by standard reasoning, the relation 

d r b'^ 



dv = - i-Bdv (B2) 

can be derived. The integral on the left-hand side is over all infinite space and thus gives the total magnetic 
energy stored in the mean magnetic field whereas that on the right-hand side is over the fiuid body only. 
Consider now the mean-field version of Ohm's law in the form 

D-j = E + VxB + 5*. (B3) 

D means the magnetic diffusivity tensor (|28|) and S* the electromotive force £ without the f3 term, 
£* = S + P ■ {V xB). The energy balance (|B2)) turns with (|B3|l into 

d r B^ 



dt 



[ ^dv = - [ {Di,jjj -]-{VxB)-]-£*)dv. (B4) 



If D is positive definite, the first term under the right integral clearly describes a decrease of the energy 
stored in the mean magnetic field. In the absence of the other two terms the field would be bound to decay. 
By the way, a part of £* with the structure S x (V x B) does not contribute to j • £*. 

We recall here that the relations (|16() have been used for the determination of (3 and so D as well as a, 
■7, S and K from the a^Ai bi^xr and ft^Ai? extracted from the numerical simulations, and that these relations 
have been derived with the particular choice bi^x^ = 0. In the magnetoconvection case D proved to be 
positive definite. Another choice for bf^x^ would lead to another D and another £*. In particular, D can 
lose its definiteness. However, the physical situation can not change by these redefinitions. £* also changes, 
such that the right-hand side of HB4|) remains its value. 

By these reasons the deviations of D from being positive definite in the geodynamo case do not seem to 
be dramatic. The growth of the magnetic energy suggested by these deviations may well be intercepted by 
induction effects described by £* . Presumably also the negative diffusivity found in the investigations by 
Brandenburg and Sokolofll (|2nn2h need not to be considered as unphysical but could also be understood in 



that sense. 



REFERENCES 

Beck, R., Brandenburg, A., Moss, D., Shukurov, A. and Sokoloff, D., Galactic magnetism: Recent devel- 
opments and perspectives. Ann. Rev. Astron. Astrophys., 1996, 34, 155-206. 

Brandenburg, A. and Sokoloff, D., Local and nonlocal magnetic diffusion and alpha-effect tensors in shear 
flow turbulence. Geophys. Astrophys. Fluid Dyn., 2002, 96, 319-344. 

Brandenburg, A. and Subramanian, K., Astrophysical magnetic fields and nonlinear dynamo theory. Phys. 
Rep., 2005a, 417, 1-209. 

Brandenburg, A. and Subramanian, K., Minimal tau approximation and simulations of the alpha effect. 
Astron. Astrophys., 2005b, 439, 835-843. 



, 2008 



13:24 Geophysical and Astrophysical Fluid Dynamics schrinner 



REFERENCES 29 

Brun, A.S., Miesch, M.S. and Toomre, J., Global-scale turbulent convection and magnetic dynamo action 

in the solar envelope. Astrophys. J., 2004, 614, 1073-1098. 
Carvalho, J.C., Dynamo effect at high order approximation. Astron. Astrophys., 1992, 261, 348-352. 
Carvalho, J.C., Diffusion effect on the dynamo action at high order approximation. Astron. Astrophys., 

1994, 283, 1046-1050. 

Christensen, U., Olson, P. and Glatzmaier, G.A., A dynamo model interpretation of geomagnetic field 

structures. Geophys. Res. Lett., 1998, 25, 1565-1568. 
Christensen, U., Olson, P. and Glatzmaier, G.A., Numerical modelling of the geodynamo: A systematic 

parameter study. Geophys. J. Int., 1999, 138, 393-409. 
Christensen, U.R., Aubert, J., Cardin, P., Dormy, E., Gibbons, S.. Glatzmaier, G.A., Grotc, E., Honkura, 

Y., Jones, C, Kono, M., Matsushima, M., Sakuraba, A., Takahashi, F., Tilgner, A., Wicht, J. and Zhang, 

K., A numerical dynamo benchmark. Phys. Earth Planet. Inter., 2001, 128, 25-34. 
Fearn, D., Hydromagnetic flow in planetary cores. Rep. Progr. Phys., 1998, 61, 175-235. 
Ferriere, K., Effect of an ensamble of explosions on the galactic dynamo. Astrophys. J., 1992, 389, 286-296. 
Ferriere, K., The full alpha-tensor due to supernova and superbubbles in the galactic disk. Astrophys. J., 

1993a, 404, 162-184. 

Ferriere, K., Magnetic diffusion due to supernova explosions and superbubbles in the galactic disk. Astro- 
phys. J., 1993b, 409, 248-261. 

Giesecke, A., Ziegler, U. and Riidiger, G., Geodynamo a-effect derived from box simulations of rotating 
magnetoconvection. Phys. Earth Planet. Inter., 2005, 152, 90-102. 

Oilman, P. A., Dynamically consistent nonlinear dynamos driven by convection in a rotating spherical shell. 
II - Dynamos with cycles and strong feedbacks. Astrophys. J. SuppL, 1983, 53, 243-268. 

Oilman, P.A. and Miller, J., Dynamically consistent nonlinear dynamos driven by convection in a rotating 
spherical shell. Astrophys. J. SuppL, 1981, 46, 211-238. 

Glatzmaier, G.A., Numerical simulations of stellar convective dynamos. I - The model and method. J. 
Comp. Phys., 1984, 55, 461-484. 

Glatzmaier, G.A., Numerical simulations of stellar convective dynamos. II - Field propagation in the 
convection zone. Astrophys. J., 1985, 291, 300-307. 

Glatzmaier, G.A. and Roberts, P.H., A three-dimensional convective dynamo solution with rotating and 
finitely conducting inner core and mantle. Phys. Earth Planet. Inter., 1995a, 91, 63-75. 

Glatzmaier, G.A. and Roberts, P.H., A three-dimensional self-consistent computer simulation of a geo- 
magnetic field reversal. Nature, 1995b, 377, 203-208. 

Kageyama, A. and Sato, T., Generation mechanism of a dipole field by a magnetohydrodynamic dynamo. 
Phys. Rev. E, 1997, 55, 4617-4626. 

Kapyla, P.J., Korpi, M.J., Ossendrijver, M. and Stix, M., Magnetoconvection and dynamo coefficients III: 
alpha-effect and magnetic pumping in the rapid rotation regime. Astron. Astrophys., 2006, in press. 

Kichatinov, L.L. and Riidiger, G., Magnetic-field advection in inhomogeneous turbulence. Astron. Astro- 
phys., 1992, 260, 494-498. 

Kowal, G., Otmianowska-Masur, K. and Hanasz, M., Dynamo coefficients in Parker instable disks with 
cosmic rays and shear. The new methods of estimation. Astron. Astrophys., 2006, 445, 915-929. 

Krause, F. and Radler, K.-H., Mean-Field Magnetohydrodynamics and Dynamo Theory, 1980 (Pergamon 
Press). 

Kuang, W. and Bloxham, J., An Earth-like numerical dynamo model. Nature, 1997, 389, 371-374. 
Kutzner, C. and Christensen, U.R., From stable dipolar towards reversing numerical dynamos. Phys. Earth 

Planet. Inter., 2002, 131, 29-45. 
McKee, S., Wall, D.P. and Wilson, S., An alternating direction implicit scheme for parabolic equations 

with mixed derivative and convective terms. J. Comput. Phys., 1996, 126, 64-76. 
Moffatt, H., Magnetic field generation in electrically conducting fluids, 1978 (Cambridge University Press). 
Nicklaus, B. and Stix, M., Corrections to first order smoothing in mean-field electrodynamics. Geophys. 

Astrophys. Fluid Dyn., 1988, 43, 149-166. 
Olson, P., Christensen, U. and Glatzmaier, G.A., Numerical modeling of the geodynamo: Mechanisms of 

field generation and equilibration. J. Geophys. Res., 1999, 104, 10 383-10 404. 



, 2008 



13:24 Geophysical and Astrophysical Fluid Dynamics schrinner 



30 REFERENCES 

Olson, P. and Glatzmaier, G.A., Magnetoconvection in a rotating spherical shell. Phys. Earth Planet. 

Inter., 1995, 92, 109-118. 
Ossendrijvcr, M., The solar dynamo. Astron. Astrophys. Rev., 2003, 11, 287-367. 

Ossendrijver, M., Stix, M. and Brandenburg, A., Magnetoconvection and dynamo coefficients: Dependence 
of the alpha effect on rotation and magnetic field. Astron. Astrophys., 2001, 376, 713-726. 

Ossendrijver, M., Stix, M., Brandenburg, A. and Riidiger, G., Magnetoconvection and dynamo coefficients. 
II. Field-direction dependent pumping of magnetic field. Astron. Astrophys., 2002, 394, 735-745. 

Parker, E.N., Hydromagnetic dynamo models. Astrophys. J., 1955, 122, 293. 

Radler, K.-H., A new turbulent dynamo. I. Monatsber. Dt. Akad. Wiss., 1969, 11, 272-279. 

Radler, K.-H., Untersuchung eines Dynamomechanismus in turbulenten Medien. Monatsber. Dt. Akad. 
Wiss., 1970, 12, 468-472. 

Radler, K.-H., Mean-field approach to spherical dynamo models. Astron. Nachr., 1980, 301, 101-129. 
Radler, K.-H., Investigations of spherical kinematic mean-field dynamo models. Astron. Nachr., 1986, 307, 
89-113. 

Radler, K.-H., Cosmic dynamos. Rev. Mod. Astron., 1995, 8, 295-321. 

Radler, K.-H., The generation of cosmic magnetic fields. In: From the Sun to the Great Attractor, 1999 
Guanajuato Lectures in Astrophysics, D. Page and J.G. Hirsch (eds.). Lecture Notes in Physics, 2000, 
101-172. 

Radler, K.-H., Apstein, E., Rheinhardt, M. and Schiiler, M., The Karlsruhe dynamo experiment - a mean- 
field approach. Studia geoph. et geod., 1998, 42, 224-231. 

Radler, K.-H., Apstein, E. and Schiiler, M., The alpha-effect in the Karlsruhe dynamo experiment. In: 
Proceedings of the International Conference "Transfer Phenomena in Magnetohydrodynamic and Elec- 
troconducting Flows" Held in Aussois, France, 1997, 9-14. 

Radler, K.-H., Kleeorin, N. and Rogachevskii, I., The mean electromotive force for MHD turbulence: 
The case of a weak mean magnetic field and slow rotation. Geophys. Astrophys. Fluid Dyn., 2003, 97, 
249-274. 

Radler, K.-H. and Rheinhardt, M., Mean-field electrodynamics: Critical analysis of various analytical 

approaches to the mean electromotive force. Geophys. Astrophys. Fluid Dyn., 2006, submitted. 
Radler, K.-H., Rheinhardt, M., Apstein, E. and Fuchs, H., On the mean-field theory of the Karlsruhe 

dynamo experiment. Nonlin. Proc. Geophys., 2002a, 9, 171-187. 
Radler, K.-H., Rheinhardt, M., Apstein, E. and Fuchs, H., On the mean-field theory of the Karlsruhe 

dynamo experiment. I. Kinematic theory. MagnetohydrodynamAcs, 2002b, 38, 41-71. 
Radler, K.-H. and Stepanov, R., The mean electromotive force due to turbulence of a conducting fiuid in 

the presence of a mean flow. Phys. Rev. E, 2006, 73, 056 311/1-15, arXiv:physics/0512120v2. 
Roberts, P.H., Kinematic dynamo models. Phil. Trans. Roy. Soc. London Ser. A, 1972, 272, 663-698. 
Rogachevskii, I. and Kleeorin, N., Electromotive force and large-scale magnetic dynamo in a turbulent 

flow with a mean shear. Phys. Rev. E, 2003, 68, 036 301/1-12. 
Riidiger, G. and Kichatinov, L.L., Alpha-effect and alpha-quenching. Astron. Astrophys., 1993, 269, 581- 

588. 

Schmitt, D., Dynamo action of magnetostrophic waves. In: Proceedings of the Fourth European Meeting 

on Solar Physics - 'The Hydromagnetics of the Sun', 1984, ESA SP-220, 223-224. 
Schmitt, D., Dynamo action of magnetostrophic waves. In: Advances in Nonlinear Dynamos, 2003 (Taylor 

and Francis), 83-122. 

Schrinner, M., Mean-field view on geodynamo models, 2005, Ph.D. thesis, Georg August University 
Gottingen. 

Schrinner, M., Radler, K.-H., Schmitt, D., Rheinhardt, M. and Christensen, U., Mean-field view on rotating 
magnetoconvection and a geodynamo model. Astron. Nachr., 2005, 326, 245 249. 

Schrinner, M., Radler, K.-H., Schmitt, D., Rheinhardt, M. and Christensen, U., Mean-field view on geo- 
dynamo models. Magnetohydrodynamics, 2006, 42, 111-122. 

Shukurov, A., On the origin of galactic magnetic fields. Astrophys. Space Sci., 2002, 281, 285-288. 

Steenbeck, M. and Krause, F., On the dynamo theory of stellar and planetary magnetic fields. I. AC 
dynamos of solar type. Astron. Nachr., 1969, 291, 49-84. 



February 5, 2008 



13:24 Geophysical and Astrophysical Fluid Dynamics schrinner 



REFERENCES 31 

Steenbeck, M., Krause, F. and Radler, K.-H., A calculation of the mean electromotive force in an electrically 
conducting fluid in turbulent motion under the influence of Coriolis forces. Zeitschr. Naturf. Teil A, 1966, 
21, 369-376. 

Tobias, S.M., The solar dynamo. Roy. Soc. Lond. Phil. Trans. Ser. A, 2002, 360, 2741-2756. 

Vainshtein, S.I. and Kichatinov, L.L., The macroscopic magnetohydrodynamics of inhomogeneously tur- 
bulent cosmic plasmas. Geophys. Astrophys. Fluid Dyn., 1983, 24, 273-298. 

Weiss, N.O. and Tobias, S.M., Physical causes of solar activity. Space Sci. Rev., 2000, 94, 99-112. 

Weiss, N., Dynamos in planets, stars and galaxies. Astrophys. Geophys., 2002, 43, 3.9-3.14. 

Wicht, J. and Olson, P., A detailed study of the polarity reversal mechanism in a numerical dynamo model. 
Geochem. Geophys. Geosyst., 2004, 5, Q03H10. 

Ziegler, U., Yorkc, H. and Kaisig, M., The role of supernovae for the galactic dynamo. Astron. Astrophys., 
1996, 305, 114-124. 



